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Abstract 

Background: Uncoupling protein one (UCPl) is a nnitochondrial inner nnembrane protein capable of uncoupling 
the electrochemical gradient fronn adenosine-5'-triphosphate (ATP) synthesis, dissipating energy as heat. UCPl plays 
a central role in nonshivering thermogenesis in the brown adipose tissue (BAT) of hibernating animals and small 
rodents. A UCPl ortholog also occurs in plants, and aside from its role in uncoupling respiration from ATP synthesis, 
thereby wasting energy, it plays a beneficial role in the plant response to several abiotic stresses, possibly by 
decreasing the production of reactive oxygen species (ROS) and regulating cellular redox homeostasis. However, 
the molecular mechanisms by which UCPl is associated with stress tolerance remain unknown. 

Results: Here, we report that the overexpression of UCPl increases mitochondrial biogenesis, increases the 
uncoupled respiration of isolated mitochondria, and decreases cellular ATP concentration. We observed that the 
overexpression of UCPl alters mitochondrial bioenergetics and modulates mitochondrial-nuclear communication, 
inducing the upregulation of hundreds of nuclear- and mitochondrial-encoded mitochondrial proteins. Electron 
microscopy analysis showed that these metabolic changes were associated with alterations in mitochondrial 
number, area and morphology. Surprisingly, UCPl overexpression also induces the upregulation of hundreds of 
stress-responsive genes, including some involved in the antioxidant defense system, such as superoxide dismutase 
(SOD), glutathione peroxidase (GPX) and glutathione-S-transferase (GST). As a consequence of the increased UCPl 
activity and increased expression of oxidative stress-responsive genes, the UCPl-overexpressing plants showed 
reduced ROS accumulation. These beneficial metabolic effects may be responsible for the better performance of 
UCPl-overexpressing lines in low pH, high salt, high osmolarity, low temperature, and oxidative stress conditions. 

Conclusions: Overexpression of UCPl in the mitochondrial inner membrane induced increased uncoupling 
respiration, decreased ROS accumulation under abiotic stresses, and diminished cellular ATP content. These events 
may have triggered the expression of mitochondrial and stress-responsive genes in a coordinated manner. Because 
these metabolic alterations did not impair plant growth and development, UCPl overexpression can potentially be 
used to create crops better adapted to abiotic stress conditions. 
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Background 

Mitochondrial uncoupling protein one (UCPl) is a nuclear- 
encoded protein located in the mitochondrial inner mem- 
brane. In the presence of fatty acids, UCPl uncouples the 
electrochemical gradient from adenosine-5 '-triphosphate 
(ATP) synthesis, dissipating energy as heat [1]. Mammalian 
UCPl has long been investigated in brown adipose tissue 
(BAT) for its role in thermogenesis and the regulation of re- 
active oxygen species (ROS) production [2-4]. 

UCPl has also been found in plants [5]. Similar to its 
mammalian orthologs, plant UCPl belongs to a multigenic 
family whose members are expressed in a time- and tissue- 
dependent manner and in response to low temperature 
[6-8]. Plant UCPs have also been shown to be involved in 
thermogenesis regulation in skunk cabbage [9] and climac- 
teric increases in respiration in fruits [10], but the wide- 
spread presence of this protein in eukaryotic organisms 
suggests that it may have other functions [8]. The overex- 
pression of Arabidopsis thaliana UCPl {AtUCPl) in to- 
bacco plants resulted in increased tolerance to oxidative 
stress [11]. In addition, the tobacco plants over expressing 
AtUCPl exhibited faster germination under control and 
stressful conditions, improved performance under drought 
and salt stresses, and increased rates of photosynthesis 
[12]. The mechanism underlying this increased stress pro- 
tection is generally associated with decreased ROS produc- 
tion [11,12], but a recent study performed on Solanum 
lycopersicum (tomato) plants overexpressing a UCP gene 
suggest a wider role for UCPs' protective mechanisms by 
altering cell redox homeostasis and antioxidant capacity 
[13]. These previous studies regarding UCP overexpression, 
along with a study of an insertional knockout of AtUCPl 
[14], suggest that UCPs may alter metabolism more glo- 
bally, modulating mitochondrial, chloroplastic and cytosolic 
metabolism. 

The role of mitochondria in energy metabolism and 
the stress response implies that this organelle commu- 
nicates with other cellular compartments. Alteration of 
mitochondrial function modulates the expression of UCPl 
and other nuclear-encoded mitochondrial proteins, includ- 
ing alternative oxidase (AOx) and type II NAD(P)H de- 
hydrogenase, through mitochondria-to-nucleus signaling 
[15,16]. Mitochondria-to-nucleus communication has been 
extensively investigated in mammalian models, in which 
a signaling process coordinates the expression of genes 
encoded by the mitochondrial and nuclear genomes 
through mechanisms known as anterograde (nucleus to 
organelle) and retrograde (organelle to nucleus) signals 
[17]. The cytosolic concentrations of Ca^^ and ATP, to- 
gether with the mitochondrial fission/fusion dynamics and 
ROS production, play central roles in the antero/retrograde 
signaling pathways and in the regulation of mitochondrial 
proliferation [18]. Both anterograde and retrograde sig- 
naling have been implicated in mitochondrial biogenesis 



[19,20], a process that is well documented in mammalian 
models [4] but poorly understood in plants [21]. In plants, 
the upregulation of nuclear-encoded mitochondrial genes 
have been seen during inflorescence development, espe- 
cially during flowering [22-24], but the molecular compo- 
nents involved in the regulation of this process remains 
elusive. Recently, it has been shown that the pentatricopep- 
tide proteins [25] and the transcription factors bZIP [26], 
WRKY [27], TCP [28] and NAC [29] may contribute to the 
retrograde regulation of mitochondrial/nucleus communi- 
cation in plants. 

In the present work, we investigated the mechanism of 
UCPl action in the stress response. We used molecular, 
cellular and genomic tools to investigate the molecular 
and cellular events resulting from the overexpression of 
AtUCPl in transgenic tobacco plants. We show that 
AtUCPl overexpression increases uncoupled respiration, 
decreases cellular ATP concentration, alters mitochondrial 
morphology, and triggers retrograde signaling, activating 
the expression of mitochondrial- and nuclear- encoded 
mitochondrial proteins. These changes are accompanied 
by a broad induction of stress-responsive genes that may 
help the A ^ZiCPi -overexpressing cells reduce the levels 
of ROS and perform better under stress conditions. The 
results are discussed in the context of the link between 
mitochondrial biogenesis and the stress response in plants. 

Results 

Overexpression oi the Arabidopsis UCPl (AtUCPl) in 
tobacco shows no phenotypic alterations independent of 
its transgene expression levels 

The AtUCPl gene, encoding the Arabidopsis ortholog of 
mammalian UCPl [6], was cloned under the control of 
the 35S CaMV promoter and transformed into tobacco 
plants [11]. Two transgenic events, P49 and P07, which 
presented intermediary and high AtUCPl expression, re- 
spectively, were chosen for this study. Although UCPl 
acts directly in the mitochondrial respiratory chain by 
uncoupling the electron transport from ATP synthesis, 
no apparent phenotypic alterations were observed in the 
P49 and P07 transgenic lines compared with the wild- 
type (WT) plants (Figure lA). This result diverges from 
previous findings showing increased shoot dry mass in 
plants overexpressing the same UCPl [12], but it is con- 
sistent with observations of tomato plants overexpress- 
ing a UCP, which did not exhibit growth stimulation 
[13]. The addition of sucrose to the growth medium may 
have altered mitochondrial metabolism and biogenesis 
[24]. In the data presented in Figure lA, we did not add 
sucrose to the nutrient solution. 

The transgenic lines P49 and P07 presented 2- to 8-fold 
increases in AtUCPl expression compared to WT plants in 
the two leaves sampled (Figure IB). Immunoblotting using 
an anti-A^UCPl polyclonal antibody, which recognizes the 
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Figure 1 Transgenic tobacco overexpressing AtUCPI exhiibits a normal plant phenotype independent of the UCP1 expression level. 

(A) Wild-type {WV) and AtUCPl transgenic lines P49 and P07 grown for 12 weeks in 50:50 soil/vermiculite mixtures fertilized with half-strength MS medium 
not supplemented with sucrose. (B) Transcript abundance of AtUCPl in LI and L3 leaves of W and transgenic lines P49 and P07. (C) Immunoblot analysis 
of LB leaves of Wand transgenic lines P49 and P07 using an ant\-AtUCPl antibody. The loading control was examined using an anti-(3-actin antibody. 
Panel (A) is representative of an experiment using 12 plants for each genotype (n = 12). Panel (B) is representative of an experiment with 6 plants for each 
genotype (n = 6). Panel (C) is representative of an experiment conducted in triplicate, ""p < 0.05 compared with the control. Error bars, mean ± s.e.m. 



tobacco and Arabidopsis UCPl, showed that the AtUCPl 
protein (approximately 32 kD) was weakly detectable in 
leaves but accumulated at high levels in P49 and P07 
transgenic lines (Figure IC). 

AtUCPl overexpression triggers mitochondria biogenesis 
through retrograde signaling 

Young and mature leaves differ in their responses to cold 
stress and photosynthetic activity [30,31]. Therefore, we 
examined whether young or more expanded leaves would 
respond differently to AtUCPl overexpression. P07, P49 
and WT plants were grown for 12 weeks, and apical leaves 
varying from youngest (LI) to oldest (L5) (Figure 2A) were 
sampled. Mesophyll protoplasts were prepared from L3, 
stained with the mitochondria-selective probe Mito- 
Tracker Red CM-H2XRos (Invitrogen, Carlsbad, CA, 
USA) and analyzed with fluorescent confocal micros- 
copy. Both P49 and P07 transgenic lines showed 1.5- to 
2-fold increases in the MitoTracker Red fluorescence 
signal compared with WT, indicating that UCPl overex- 
pression led to an increase in the mitochondrial number 
and/or volume (Figure 2B, C, and D). The Arabidopsis 
mitochondrial DNA (mtDNA) encodes 57 proteins 
involved in mitochondrial metabolism and mtDNA 
maintenance [32]. In mammals, mtDNA replication is 
strongly associated with cellular demand for ATP [33]. 
To investigate whether the increased MitoTracker Red 



fluorescence signal observed in the P49 and P07 trans- 
genic lines is associated with alteration in mtDNA, we 
used the MATR maturase gene to quantify its mtDNA 
content relative to the WT line. qRT-PCR quantification 
of MATR in total DNA extracted from LI to L5 leaves re- 
vealed an estimated mtDNA increase of 2- to 5-fold in the 
P49 and P07 transgenic lines, respectively, compared with 
that of WT plants (Figure 2E). Because mitochondrial 
biogenesis requires the coordinated expression of both 
nuclear and mitochondrial genomes [17,24], we exam- 
ined whether AtUCPl overexpression results in the up- 
regulation of other nuclear- and mitochondrial-encoded 
mitochondrial proteins. q-PCR of total RNA extracted 
from LI to L5 leaves from P49, P07 and WT plants was 
used to quantify relative mRNA content from genes encod- 
ing a selected set of nuclear- and mitochondrial-encoded 
mitochondrial proteins. Independent of leaf age, the P49 
and P07 transgenic lines presented 2- to 5-fold increases 
in the transcript abundance of NADH dehydrogenase 
(NADH-DeH) and NADH iron-sulfur proteins 2 and 7 
(NADH IS2-7), which are components of mitochondrial 
respiratory chain complex I [34] (Figure 2F). The increases 
in both nuclear-encoded (NADH-DeH and NADH-IS7) 
and mitochondrial-encoded (NADH-IS2) components of 
complex I indicate that mitochondrial/nuclear coordination 
is necessary for bioenergetic adaptation, mincreased tran- 
script abundance in P49 and P07 transgenic lines was also 
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Figure 2 Overexpression of AtUCPl in tobacco induces mitochondrial biogenesis through retrograde signaling. (A) Five leaves, from 



youngest (LI) to oldest (L5), were sampled from WT and transgenic lines P49 and P07. (B and C) Mesophyll protoplasts isolated from WT (B) and 
P07 (C) L3 leaves stained with MitoTracker Red. The panels are representative of 50 protoplasts evaluated for each genotype. Scale bars = 5 pm. 
(D) Quantification of the MitoTracker Red signal in a population of protoplasts isolated from the L3 leaves of WT and transgenic lines P49 and 
P07 (n = 50). (E) Relative mtDNA content of W and transgenic lines P49 and P07 in leaves LI to L5. (F) Relative mRNA content of genes involved 
in targeting proteins to mitochondria (GRPe chaperone) and energy metabolism (NADH-DeH, NADH-IS2, and NADH-IS7) in LI leaves of WT and 
transgenic lines P49 and P07. (G) Immunoblot analysis of nuclear- (IDH, VDAC and AOX) and mitochondrial-encoded (COXII) proteins involved in 
energy metabolism, transport, and stress response in the LI leaves of WT and transgenic line P07. The data presented in panels E and F are 
representative of four biological replicates (n = 4). ""p < 0.05 and """"p < 0.1 compared with WT. Error bars, mean ± s.e.m. 



observed for the molecular chaperone GRPe and for the 
adenine nucleotide exchange factor of heat shock protein- 
70 (HSP70) that is involved in mitochondrial targeting 
of proteins [35] (Figure 2F and Additional file 1: Figure 
SIA-D). We also observed that the P07 transgenic line 
showed 2- to 5-fold increases in the protein levels of 
cytochrome oxidase subunit II (COXII), isocitrate de- 
hydrogenase (IDH), voltage-dependent anion channel- 1 
(VDACl), and alternative oxidase (AOX), which are pro- 
teins associated with both oxidative stress [36] and mito- 
chondrial metabolism [37-39] (Figure 2G). These data 
indicate that AtUCPl overexpression triggered retrograde 
signaling, leading to increased transcription and trans- 
lation of mitochondrial proteins encoded by both the 
nuclear and mitochondrial genomes, culminating in in- 
creased mitochondrial biogenesis. 

AtUCPl overexpression increases uncoupling respiratory 
capacity and decreases cellular ATP content 

Although AtUCPl overexpression protects tobacco plants 
from oxidative, osmotic and drought stress, no previous 
data indicate whether UCPl does, in fact, uncouple ATP 
synthesis in these plants [11]. However, previous work has 
demonstrated that mitochondria isolated from Solanum 
tuberosum (potato) plants overexpressing a UCP show 



increased uncoupled respiration [40]. Here, we demon- 
strate that crude mitochondria isolated from L3 leaves of 
WT and P07 plants retained their respiratory proper- 
ties: ADP and FCCP can efficiently increase oxygen con- 
sumption when compared to basal respiration (Figure 3A). 
Additionally, transgenic plants showed decreased adeno- 
sine 5 '-diphosphate (ADP)-dependent respiration and 
increased uncoupling respiration in comparison to WT 
plants (Figure 3B). There was no significant difference 
in the oxygen consumption in the presence of carbonyl 
cyanide-4-trifluoromethoxy phenylhydrazone (FCCP) 
between P07 and WT mitochondria (Figure 3B). The in- 
creased uncoupling respiration of P07 compared with WT 
mitochondria is evident in a comparison of the leak ratios 
(uncoupled/maximum oxygen consumption) in P07 and 
WT (Figure 3B). The increased uncoupling respiration in 
the P07 transgenic line resulted in 20-35% decreased cel- 
lular ATP concentration in the L3 leaves compared with 
the WT (Figure 3C). 

Altered mitochondrial metabolism due to Increased UCP1 
activity leads to changes In mitochondrial number, 
volume and morphology 

Mitochondrial fission and fusion dynamics, as well as 
mitochondrial morphology, have been shown to be directly 
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Figure 3 Overexpression of AtUCPI increases mitochondrial uncoupled respiration, resulting in decreased cellular ATP content. 

(A) Time-course of respiration measurements on isolated mitocliondria from WT and P07 plants. Malate (10 mM) and glutamate (5 mM) were 
added to the respiration medium to provide substrates for the electron transport chain. ADP (800 pM) was used to stimulate ADP-dependent 
respiration, and the ATP synthase inhibitor oligomycin (5 pM) was used to measure oxygen consumption due to uncoupled respiration. Maximum 
oxygen consumption was measured in the presence of FCCP (2 pM). (B) Average oxygen consumption ofWT and P07 isolated mitochondria in 
the presence of ADP, oligomycin and FCCP. The leak ratio was calculated as the ratio between oxygen consumption in the presence of oligomycin 
and that in the presence of FCCP. (C) Decreased ATP levels in the L3, L4 and L5 leaves of the P07 transgenic line compared with WT. The data are 
representative of six biological replicates (n = 6). ""p < 0.05 and **p < 0.1 compared with WT. Error bars, mean ± s.e.m. 



affected by modulations in the energy demand and nutri- 
ent supply in human cells [41]. The fine regulation of these 
parameters is linked to the adaptation of mitochondrial 
architecture to metabolic demand. To address whether 
AtUCPl overexpression affects mitochondrial morphology, 
we performed transmission electron microscopy (TEM) of 
mesophyll cells of L2 leaves from P07 and WT plants. Ap- 
proximately 20 micrograph fields for each genotype were 
analyzed for mitochondrial number and area using ImageJ 
[42]. Images from P07 leaves revealed increases of 1.6-fold 
in mitochondrial number and 1.4- fold in mitochondrial 
volume, resulting in an overall augmentation of 2.2-fold in 
the total mitochondrial area (Figure 4A, B, and C). Surpris- 
ingly, and perhaps for the first time, we observed a novel 
mitochondrial morphology in plant cells: donut-shaped 
mitochondria (Figure 4D). This morphology has been 
observed in human cells during treatment with FCCP or 
subjection to hypoxic stress [43]. These treatments decrease 
oxygen availability for ATP production. This unusual mor- 
phology has also been observed during reoxygenation and 



FCCP washout, when cells reoxygenate and ATP levels 
are partially restored. When ATP levels were completely 
restored, the donut-shaped mitochondria fragmented, 
generating new mitochondria. These findings suggest 
that AtUCPl overexpression triggers a retrograde signal 
that promotes mitochondrial biogenesis. The alterations 
in both mitochondrial architecture and morphology may 
be due to metabolic stress imposed by the higher activity 
of the electron transport chain, thus increasing oxygen 
consumption, together with the reduced ATP synthesis 
caused by the AtUCPl overexpression. The donut-shaped 
mitochondrial morphology might be associated with the 
fusion/fission dynamics but may also be important to alle- 
viate the tension imposed by the increased mitochondrial 
volume [43]. 

AtUCPl overexpression induces broad transcriptomic 
activation of mitochondrial and stress-responsive genes 

It was unexpected that the overexpression of a single pro- 
tein would result in mitochondrial biogenesis because the 
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Figure 4 Overexpression of AtU CP 1 increased mitochondrial biogenesis through donut formation/fragmentation. (A) Representative TEM 
micrograph of WT L2 leaves. (B) Representative TEM micrograph of P07 L2 leaves. (C) Quantification of the mitochondrial number, area, and total 
area occupied by mitochondria in TEM micrograph field images of Wand P07 plants. (D) Close-up of the donut-shaped mitochondria at different 
developmental stages in P07 transgenic cells. Image analysis was performed on 20 field images (n = 20) for each genotype, ""p < 0.05 compared with 
WT. Error bars, mean ± s.e.m. Scale bars, 1 pm. 



mitochondrial proteome is composed of over 2000 pro- 
teins [24]. Therefore, the induction of mitochondrial bio- 
genesis would require the expression of an array of genes 
encoding these proteins. To examine the extent of retro- 
grade signaling activation due to AtUCPl overexpression, 
we conducted a global transcriptomic analysis of L3 leaves 
from P07 and WT plants. The transcripts were annotated 
using UniProt [44] and categorized into the Clustering of 
Orthologous Groups (COG) database [45]. We used TAIR 
subcellular prediction of proteins [46] to identify a total of 
705 nuclear-encoded mitochondrial transcripts that have 
significantly higher RPKM values in the P07 transgenic 
line and 165 that have higher expression levels in WT 
plants (Figure 5A, red dots). In addition, we identified 12 
mitochondrial-encoded proteins, all of which were more 
than 2-fold (significantly) upregulated in the P07 trans- 
genic line (Figure 5 A, blue dots). The most representative 
mitochondrial upregulated proteins include those associ- 
ated with energy production, pentatricopeptide proteins, 
proteins associated with mitochondrial transcription and 
translation, and a large number of unknown proteins 



(Figure 5B and Additional file 2: Table SI). Proteins with 
unknown subcellular prediction were submitted to further 
analysis using TargetP [47]. This analysis allowed us to 
identify another 43 proteins predicted to be targeted to 
mitochondria with fold-changes higher than 2.0 in P07 
plants (Additional file 3: Table S2). Among several proteins 
with unknown function, we identified possible orthologs of 
one of the mitochondrial calcium uniporters (MCUs) and 
a mitochondrial sodium/calcium exchanger (MCX). One 
of these proteins (Solyc04g0799 10.2.1 in Additional file 3: 
Table S2), which is 2.5-fold upregulated in P07, con- 
tains a coiled-coil domain on the C-terminal region, 
which is found among the MCUs, and has a strong se- 
quence similarity with the Arabidopsis thaliana (A^MCU), 
Mus musculus (MmMCU) and Homo sapiens (HsMCU) 
mitochondrial calcium uniporters. The other protein (Solyc 
07g042000.2.1 in Table 2) contains a sodium/calcium- 
exchanging domain whose sequence is predicted to be 
targeted to mitochondria. Among the upregulated pro- 
teins in the P07 transgenic line, we also identified several 
proteins associated with calcium signaling, mitochondrial 



Barreto et a I. BMC Plant Biology 2014, 14:144 
http://www.bionnedcentral.conn/1471-2229/14/144 



Page 7 of 1 5 



10 



^ 4 
Q_ ^ 

CM 

o 

0 

o 

CL 

-2 
-4 
-6 



— Median 




• MEP 


. • ' ''ww* •• 


• NEMP 
TP 




























n *•••., 





CL 

q: 

CM 

O 



o 
Q_ 



-2 0 2 4 6 
WT (Log2 RPKM) 



8 10 



10 
8 
6 
4 
2 
0 
-2 
-4 
-6 



— ivieaian 






• SRP 

r JP 








• *•* t*Ji 


















'****. 
















/ n *.••*. 







B 



-2 0 2 4 6 
WT (Log2 RPKM) 



8 10 



o 

o 



Energy production 
Pentatricopeptide proteins 
Transcription and translation 
Lipid metabolism 
Aminoacid metabolism 
Auxin-induced proteins 
Other function 
Unkwnon function 



60i 
45 
30 
15 

5t| 

4 
3 

2-1 
1 

0^ 

Figure 5 Transcriptome analysis reveals a global upregulation of mitochondrial and stress-responsive genes in transgenic P07 compared 
with WT. (A) Scatterplot of log2 RPKM showing nuclear- (Red circle symbol) and mitochondrial (Blue circle symbol)-encoded mitochondrial proteins 
that were upregulated in the P07 AtUCPl transgenic line compared with WT. The green dots represent the total differentially expressed genes. 
(B) Examples of genes associated with nuclear and mitochondrial-encoded mitochondrial proteins that were upregulated in the P07 transgenic line 
compared with WT. (C) Scatterplot of log2 RPKM showing stress-responsive (Red circle symbol) genes that were upregulated in AtUCPl transgenic lines 
compared with WT. The green dots represent the total differentially expressed genes. (D) Examples of genes associated with abiotic stress responses 
that were upregulated in the P07 transgenic line compared with WT. Normalization was performed by applying log2 to the RPKM values of W and 
the P07 transgenic line. All the represented transcripts showed significantly different expression between the genotypes (p < 0.05). The RPKM means 
were calculated using four biological replicates (n = 4) of L3 leaves from each genotype. 





import machinery and mitochondrial fission and fusion 
processes, along with transcripts related to lipid metabol- 
ism, which may act to support UCPl activity and cell me- 
tabolism (Additional file 4: Table S3). Several transcription 
factors that may play a role in mitochondrial retrograde 
signaling were also identified (Additional file 4: Table S3). 
Surprisingly, we also identified 1071 genes responsive to 
abiotic stress, 72% (770) of which have increased expression 
in the P07 transgenic line compared with WT (Figure 5C, 
red dots). Among these stress-associated proteins, we iden- 
tified large numbers associated with heat, osmotic, cold, 
drought, salt, oxidative and cadmium stress (Figure 5D). 
We also observed 1.4- to 3.4-fold increases in the transcript 



levels of genes involved in antioxidant mechanisms, includ- 
ing SOD, GPX and GSTs, among the upregulated stress- 
responsive genes (Additional file 5: Table S4). 

UCPl transgenic line performed better under 
abiotic stresses 

In addition to the previous observation of increased tol- 
erance to oxidative, salt, and drought stresses [11,12], we 
observed that the P07 UCPl overexpressor performed 
better under low and high pH, hyperosmotic stress, and 
oxidative stress (Figure 6A). Additionally, in the presence of 
free fatty acids, the P07 transgenic line also exhibited im- 
proved performance under low temperature, a result that 
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Figure 6 >AfL/CP7-overexpressing plants performed better under abiotic stress, most lilcely by decreasing ROS accumulation. 

(A) Two-week-old WT and P07 plantlets were plated on half-strength MS medium supplemented with 80 mM mannitol (hyperosmotic stress) or 



10 mM H2O2 (oxidative stress), pH 7.5 and pH 4.5 for 22 days under a day/light regime of 16 h at 100 m 



and Shin darkness. For the 



low-temperature treatment, the seeds were germinated and grown for 30 days on half-strength MS medium supplemented with 50 |jM:10 |jM 
linoleic acid:bovine serum albumin. Plant growth assays were performed using four biological replicates (n = 4) on plates containing 5 plants for 
each treatment. *p < 0.05 and """"p < 0.1 compared with WT. Error bars, mean ± s.e.m. Scale bars, 200 |jm. (B) Leaf discs (0.8 cm) or (C) root tips 
isolated from WT and the P07 transgenic line were incubated in half-strength MS medium supplemented with 100 mM NaCI, lOOmM mannitol, 
or 10 mM H2O2 for 20 h. Control treatments were performed using unsupplemented half-strength MS medium. Tissues were stained with MitoSox 
red and visualized with fluorescent confocal microscopy. 



support the thermogenic role of plant UCPl (Figure 6A), 
Because UCPl activity in mouse brown adipose tissue, po- 
tato tuber and tobacco leaves can reduce superoxide accu- 
mulation [12,40,48], we examined whether this process also 
occurs in the P07 transgenic line but used an approach that 
detects ROS accumulation only within mitochondria. When 
subjected to salt, oxidative and hyperosmotic stress, mito- 
chondrial ROS accumulation was decreased in the P07 



transgenic line compared with WT plants in both leaf discs 
and roots (Figure 6B and C). 

Discussion 

The mitochondrial respiratory chain is composed of four 
macromolecular complexes through which free energy is 
conserved by coupling electron transport to the forma- 
tion of a proton motive force by complexes I, III, and IV, 
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which is then dissipated by FlFO-ATPase (complex V) 
for ATP synthesis [49]. In the BAT, electron transport is 
uncoupled from ATP synthesis by UCPl for heat produc- 
tion [4]. However, BAT is a specialized tissue dedicated to 
thermogenesis regulation whose main characteristic is a 
dramatic increase in mitochondria number and UCPl 
expression [50]. In other cell types, however, decreased 
ATP concentration may compromise cell metabolism [51]; 
therefore, cells somehow sense the cytoplasmic ATP con- 
centration as a signal to change gene expression to quickly 
readjust the mitochondrial architecture and cell metabol- 
ism [18,52]. In this work, we show that the overexpression 
of UCPl increases uncoupled respiration and decreases 
cellular ATP concentration. There is evidence demonstrat- 
ing that mitochondria isolated from potato plants overex- 
pressing UCP exhibit increased uncoupling respiration 
[40], but the present study represents the first time that 
the effect of UCP on intracellular ATP level has been 
quantified in plants. It is possible that the enhanced un- 
coupled respiration and its consequent negative effect on 
the ATP production trigger a signal to compensate for the 
energy metabolism and maintain ATP at safe levels. Con- 
sistent with this view is the remarkable upregulation of 
electron transport chain components, including NADH- 
DeH, NADH-IS2-7, and COX, observed in the P07 over- 
expressor. The coordinated upregulation of the electron 
transport chain components encoded by the nuclear 
and mitochondrial genomes occurs in tissues with high 
metabolic demands and increased mitochondrial bio- 
genesis [21,50,53], such as flowers [22-24,54] and ger- 
minating seeds [55]. 

Mitochondrial biogenesis was observed using multiple 
approaches, including MitoTracker RED fluorescence 
imaging, quantitation of mtDNA and TEM. In all cases, 
these methods allowed us to conclude that the overex- 
pression of AtUCPl resulted in increased mitochondrial 
number and volume. To assess whether the upregulation 
of mitochondrial genes was restricted to the electron 
transport chain components, we performed global tran- 
scriptome sequencing and observed the upregulation of 
hundreds of genes encoding mitochondrial proteins. All 
of the identified transcripts encoded on the mtDNA were 
upregulated in transgenic plants, which is consistent with 
the increased mtDNA content. For instance, the ATP syn- 
thase gene encoded by the mtDNA was found to be up- 
regulated by 8.1 -fold, and two of its nuclear subunits were 
increased by 2.7- and 3.6-fold. This upregulation was not 
restricted to the components involved in the oxidative 
phosphorylation but also included genes that participate 
in the tricarboxylic acid cycle (TCA). The key control 
point of the TCA is the enzyme IDH [56], which showed 
increased protein levels in P07 plants (Figure 2G). In fact, 
experiments performed with mitochondria isolated from 
potato plants demonstrated that UCPl activity facilitates 



higher TCA flux by decreasing membrane potential and, 
consequently, ROS production [40] . 

Transcriptional regulation is the main regulatory factor 
in the expression of nuclear-encoded genes associated 
with respiration [21], but the vast majority of proteins 
targeted to the mitochondria are transported post- 
translationally by tightly regulated machinery [57]. Among 
these proteins, HSP70 is involved in the transport from the 
cytoplasm to the mitochondrial outer membrane, whereas 
the mitochondrial TOM and TIM translocases are re- 
sponsible for transport inside the organelle [58]. Inter- 
estingly, all of these proteins were upregulated in the P07 
AtUCPl overexpression line. The increased expression of 
components of the mitochondrial energy production ma- 
chinery, along with the increased mitochondrial number of 
transgenic cells, may alter cellular energy homeostasis. 
Therefore, the mitochondria may consume more substrates 
in an attempt to boost ATP synthesis; for that purpose, 
these substrates must be supplied from other sources. This 
need might explain the increased photosynthetic rate 
exhibited by these plants under control conditions [12]. 
Transcriptome sequencing provides an important tool 
to investigate this possibility. We did not observe differ- 
ential expression patterns between WT and P07 plants 
regarding chloroplast- targeted proteins (data not shown), 
but further analysis of the transcriptome might provide 
new insights into this issue. 

The adaptation to changes in energy supply and de- 
mand influences not only mitochondrial gene expression 
but also mitochondrial architecture and morphology [41]. 
Alterations in mitochondrial morphology during starva- 
tion have already been observed in Arabidopsis thaliana 
cells cultured without a carbon source [24]. In human 
cells, treatment with the chemical uncoupler FCCP causes 
an increase in oxygen consumption and a decrease in 
ATP generation, provoking mitochondrial fragmentation 
and the appearance of donut-shaped mitochondria [43]. 
This fragmentation may occur to support an increase in 
respiration in an attempt to elevate ATP levels [41], 
whereas the formation of the donut-shaped morphology 
may be a component of a protective mechanism that helps 
to preserve the organelles from damage under conditions 
of metabolic stress [41,43]. When FCCP is washed out of 
these cells, the donut morphology persists until the ATP 
levels are completely restored, at which point the donuts 
fragment and new mitochondria arise. The mechanisms 
underlying the donut formation in UCPl -over expressing 
cells might be explained by an analogy with human cells. 
Increased UCPl activity decreases the membrane poten- 
tial that is used for ATP synthesis, provoking an increase 
in mitochondrial number to balance the ATP levels in 
transgenic plants. Mitochondria are responsible for at 
least 90% of cell oxygen consumption [59]. We therefore 
expect that an increase in respiratory chain machinery 
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and mitochondrial number, together with the constant 
leak caused by UCPl, would result in increased oxygen 
consumption per molecule of ATP generated. Thus, this 
metabolic stress caused by UCPl overexpression appears 
sufficient to allow the formation of donut mitochondria in 
plant cells. 

Communication between mitochondria and the nu- 
cleus is required for the adaptations that we observed in 
P07 transgenic plants. Although several important regu- 
lators of chloroplastic retrograde signaling have been ex- 
tensively studied [60], relatively little is known about the 
regulators of mitochondrial retrograde signaling in plants. 
Several recently published studies suggest roles for the 
WRKY, bZIP, TCP, and NAC transcription factor families 
in activating the expression of nuclear- encoded mitochon- 
drial proteins in plants [26-29]. These transcription factors 
recognize elements that are overrepresented in the regula- 
tory region of mitochondrial genes and that are involved 
in several biotic and abiotic stress responses [21]. The 
colocalization of proteins in mitochondria and the nucleus 
appears to be an important topic that is very little ex- 
plored in plants. For example, an Arabidopsis thaliana 
pentatricopeptide protein (PNMl) [25] was identified in 
both mitochondria and nucleus that may act to coordinate 
the expression between the two genomes. The increased 
expression of transcription factors that belong to the 
WRKY and TCP families, along with the upregulation of a 
large number of pentatricopeptide proteins in transgenic 
plants, presents the opportunity to further explore new 
candidates involved in mitochondrial retrograde signal- 
ing. The expression of most transcriptional regulators 
is further linked to cellular signals and environmental 
cues [21]. The most studied cellular signals involved in 
mitochondria-nucleus signaling are ATP, calcium and 
ROS [18]. Although the role of mitochondria in calcium 
metabolism has been extensively studied in mammals, lit- 
tle is known about calcium signaling by plant mitochon- 
dria [61]. Mitochondria in mammals are able to sequester 
cytoplasmic Ca"^^, functioning as a transient calcium store 
for protective mechanisms [18,62]. Additionally, the Ca"^^ 
levels inside the organelle are positively correlated with in- 
creased ATP production [63], although Ca^^ overaccumu- 
lation is linked to the induction of apoptosis in both 
mammals and plants [61]. Consistent with this relation, 
it is very interesting that we found two genes upregu- 
lated in transgenic plants that are possible orthologs, 
an MCU and an MCX, along with several genes in- 
volved in calcium signaling (Additional file 3: Tables S2 
and Additional file 4: Table S3). The MCUs and MCXs 
have been extensively studied in animal cells, and their 
expression is tightly linked with mitochondrial calcium 
transport [64,65]. 

The main focus of previous studies regarding UCPl over- 
expression in plants has been ROS production [11-13,40]. 



Both potato and tobacco transgenic plants presented 
decreased superoxide production under abiotic stresses 
[11,12,40]. Additionally, tomato plants overexpressing 
UCPl demonstrate not only an alteration on ROS pro- 
duction but also a possible role of UCPl in regulating 
the cellular redox homeostasis [13]. The plants analyzed 
in the present study showed increased expression not 
only of the antioxidant defense system, including GPX, 
GSTs and SOD, but also of an array of stress-responsive 
genes. This effect may help explain the positive impact 
of AtUCPl overexpression on the tolerance to multiple 
abiotic stresses. Consistent with previous data, we show 
in this work that AtUCPl overexpression protects plants 
from a number of abiotic stresses, including high pH, low 
pH and cold. The fact that transgenic tobacco performs 
better at low temperatures when the medium is supple- 
mented with fatty acids reinforces the thermogenic role of 
UCPl in plants. The importance of ROS in cell signaling 
is well known [66]; thus, the results regarding decreased 
mitochondrial ROS production under oxidative, osmotic 
and salt stresses also reinforce the importance of this mol- 
ecule in UCPl -mediated retrograde signaling. 

Together, these data indicate that A^ZiCPi -overexpressing 
lines suffer metabolic stress caused by the increased un- 
coupling respiration. Chronic overexpression of UCPl can 
induce mitochondrial biogenesis in mammals [67]. There- 
fore, we propose that in the transgenic plants, the lower 
ATP level acts as a key element in the retrograde signaling, 
promoting a broad increase in the expression of both mito- 
chondrial genes and stress-related genes. This pattern is 
reflected in the decreased ROS production and the better 
performance of these plants under various stresses. These 
findings advance our understanding of stress -tolerance 
mechanisms, mitochondrial biogenesis and bioenergetic 
adaptation in plants and therefore might assist in the 
implementation of biotechnological tools for the develop- 
ment of abiotic stress -tolerant plants. The overall regulation 
of stress-responsive genes indicates a link between UCPl- 
activated mitochondrial biogenesis and the increased stress 
response exhibited by the transgenic plants. 

Conclusions 

In this work, we present strong evidence that the over- 
expression of UCPl protein affects the mitochondrial 
dynamics at both structural and metabolic levels, leading 
to an increased mitochondria number and volume per 
ceU. It could be argued that these changes are a conse- 
quence of the overpopulation of UCPl protein at the 
inner mitochondrial membrane, thereby leading to abnor- 
mal mitochondrial structure and function. We interpret the 
results presented in this manuscript as a beneficial effect of 
the overexpression of UCPl that triggers a retrograde 
signaling process, signaling the nuclear and mitochondrial 
genomes to increase the production of mitochondrial 
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proteins. Consequently, more mitochondrial membrane 
is produced. The increase in UCPl, both because of the 
overproduction of the protein and the increased mitochon- 
drial number and volume, leads to increased uncoupled 
respiration, decreasing the cellular ATP concentration. This 
change may trigger metabolic stress and enhances a strong 
stress response. Because it does not affect plant growth and 
development, this mechanism might be used to create 
crops better adapted to abiotic stress conditions. 

Methods 

Plant materials and growth conditions 

Tobacco {Nicotiana tabacum) ecotype SRI overexpressing 
AtUCPl was obtained as previously described [11]. In this 
work, we used two independent transgenic lines, P49 and 
P07, which express AtUCPl at intermediate and high levels, 
respectively. Transgenic and wild-type (WT) plants were 
grown under a light/dark cycle of 16 h at 100 (lE m"^ s"^ 
and 8 h in darkness for 12 weeks on 50:50 soil/vermiculite 
mixtures at 24°C. The plants were fertilized weekly with 
half-strength MS medium not supplemented with sucrose. 
For the cellular and molecular analysis, samples of the five 
youngest leaves (LI to L5) of each genotype were collected, 
frozen in liquid nitrogen, and stored at -80°C. For the other 
experiments, seeds were grown directly in experimental 
medium, as described in the figure legends. 

qRT-PCR 

For mitochondrial DNA quantification, total DNA was 
isolated from WT, P49, and P07 plants using the plant 
DNAzol reagent (Invitrogen, Carlsbad, CA, USA) accord- 
ing to the manufacturer s protocol. For the mRNA expres- 
sion analysis, total RNA was isolated from leaf tissue using 
the NucleoSpin Mini RNA/Protein Kit (Macherey-Nagel, 
Duren, Germany). cDNA was synthesized using the Rever- 
taid First Strand cDNA Synthesis Kit (Fermentas, Vilnius, 
Lithuania) according to the manufacturer s protocol. qRT- 
PCR was performed using an ABI PRISM 7500 (Applied 
Biosystems, Foster City, CA, USA) with SYBR Green 
(AppUed Biosystems). The reactions were performed at 
least in triplicate for four biological replicates. The mRNA 
expression and mtDNA content analysis were normalized 
using ACTINl as an internal reference. Values are pre- 
sented as the mtDNA content or transcript abundance of 
P49 and P07 relative to WT. Students t-test was per- 
formed to determine significance (p < 0.05). Primers based 
on the MATR gene were used for mitochondrial DNA 
quantification [68]. 

The following primers were designed for amplification: 

ACTINl 

FW: 5 '-ACTGTCCACGAGGTCCGG-3' 
RV: 5'-TGTCGGATCTTGCGCGGC-3' 
AtUCPl 

FW: 5'-TTGAGCAAGAAAATTCTTGCTG-3' 



RV: 5'-AGGCGGAAGGAAAATTAGC-3' 
NADH-DeH 

FW: 5'-GGCTGAAGCGCGAGAAGAC-3' 
RV: 5'-CAGGGCAGGCTTCTTGGC-3' 
NADH-IS2 

FW: 5'-GTGAAGTGGCGTGGCAAAC-3' 
RV: 5'-TTGTGGATCGCGAAGGGAG-3' 
NADH-IS7 

FW: 5'-CGCAGTTGCGAAGCGAAAC-3' 
RV: 5'-GGGCCATCGCGACAGAGA-3' 
GRPe chaperone 

FW: 5'-AAACCTTGGCTTGTGACCCA-3' 
RV: 5'-TCATTCGGCCAGCTGAAGTT-3' 

Immunoblot analysis 

Total soluble proteins were extracted from WT, P49 and 
P07 plants by grinding frozen leaf tissue in 5x protein 
sample buffer (40 mM Tris, pH 7.4; 30 mM NaCl; 
10 mM p-mercaptoethanol; 0.1% Triton X-100; 5 mM 
benzamidine). The proteins were separated on 12% SDS- 
PAGE gels, transferred to nylon membranes, and blotted 
using the following polyclonal antibodies raised against 
Arabidopsis proteins (Agrisera, Vannas, Sweden): COXII, 
VDAC (At3g01280), AOxl/2 (At2g22370 and At5g64210), 
IDH (At4g35260 and At2g27130), and ACTIN (At2g37620, 
At3gl8780, At3g53750, At5g59370, At2g42100, At5g09810, 
Atlg49240, At3gl2110, and At3G46520). The bands were 
detected using an ImageQuant LAS500 (GE Healthcare, 
Little Chalfont, UK) with the SuperSignal West Pico 
chemiluminescent substrate (Thermo-Scientific, Waltham, 
MA, USA). 

Mitochondrial quantification in leaf mesophyll protoplasts 
by MitoTracker Red staining 

Tobacco mesophyll protoplasts were isolated from WT, 
P49, and P07 plants as previously described [69], with 
minor modifications. Leaf strips (0.1 cm in length) were 
vacuum- infiltrated and incubated for 3 h in an enzyme 
solution (20 mM MOPS, pH 5.7; 0.4 M d-mannitol; 20 mM 
KCl; 10 mM CaC^; 0.1% BSA; 1.5% cellulase; 0.4% macero- 
zyme). After digestion, the suspension was filtered through 
75 (im nylon mesh and centrifuged at 100 x g for 2 min. 
The protoplasts were resuspended in W5 buffer (2 mM 
MOPS, pH 5.7; 154 mM NaCl; 125 mM CaC^; 5 mM 
KCl), incubated on ice for 30 min, and subsequently resus- 
pended in MMG solution (4 mM MOPS, pH 5.7; 0.4 M 
mannitol; 15 mM MgCl2). The protoplasts were then 
stained with 250 nM MitoTracker Red CM-H2XRos 
(Invitrogen) for 30 min in the dark and washed in 
MMG solution before viewing. MitoTracker Red CM- 
H2XRos fluorescence was detected at Ex/Em of 543/ 
589 nm and chlorophyll auto fluorescence was detected 
at Ex/Em 543/645 nm using a confocal microscope 
(LSM780-NLO; Zeiss, Oberkochen, Germany) equipped 
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with a 40x oil-immersion objective. The image analysis 
was performed for 40 protoplasts per sample using ImageJ 
1.44p [42]. MitoTracker Red CM-H2XRos fluorescence 
intensity was normalized using the chlorophyll autofluo- 
rescence intensity. 

Mitochondrial superoxide detection by MitoSox staining 

Leaf discs (0.8 cm in diameter) isolated from L3 leaves were 
treated with half-strength MS medium (control) or half- 
strength MS medium supplemented with 100 mM NaCl 
(salt stress) or 10 mM H2O2 (oxidative stress) for 20 h 
under constant light and agitation. For root analysis, plants 
grown for 3 weeks on half-strength MS medium were 
transferred to petri dishes containing fresh half-strength 
MS medium not supplemented (control) or supplemented 
with 100 mM NaCl (salt stress) or 100 mM mannitol 
(hyperosmotic stress) and kept under constant light for 
20 h. The leaf discs and roots were stained with 5 [iM 
MitoSox-Red (Invitrogen) for 25 min and washed with half- 
strength MS medium before viewing. Fluorescence was 
detected at Ex/Em of 510/580 nm using a confocal micro- 
scope (LSM780-NLO; Zeiss) equipped with a lOx objective. 

Transmission electron microscopy 

Leaf strips (1-2 mm in length) from L2 leaves of WT 
and P07 plants were submerged and vacuum infiltrated 
for at least 2 h in fixative solution (2.5% glutaraldehyde, 
2.5% paraformaldehyde in 0.1 M sodium cacodylate buffer, 
pH 7.4). The samples were washed in 0.1 M cacodylate buf- 
fer and postfixed in a solution containing 1% osmium tet- 
roxide (OSO4) and 1.5% potassium ferrocyanide (KFeCNe) 
for 1 h. The samples were then washed 3 times in water 
and incubated in 1% uranyl acetate in maleate buffer for 
1 h; this was followed by 3 washes in maleate buffer and 
subsequent dehydration in a graded ethanol series (50%, 
70%, and 90%, 10 min each; 100%, 2 x 10 min). The sam- 
ples were then placed in propylene oxide for 1 h and infil- 
trated in a 1:1 mixture of propylene oxide and TAAB 812 
Resin mixture (Marivac Canada, St. Laurent, Canada). The 
samples were embedded in TAAB 812 Resin mixture and 
polymerized at 60°C for 48 h. Ultrathin sections of approxi- 
mately 90 nm were cut using a Reichert Ultracut-S micro- 
tome, placed onto copper grids, stained with lead citrate, 
and examined using a Tecnai G2 Spirit BioTWIN transmis- 
sion electron microscope. The images were recorded with 
an AMT 2 k CCD camera. The mitochondrial number 
and area were quantified on 20 field images representing 
each genotype at 4800 x direct magnification. The image 
analysis was performed using ImageJ 1.44p [42]. 

ATP measurement 

ATP assays were performed on L4 and L5 leaves isolated 
from four biological replicates of WT and P07 plants. ATP 
extraction was performed as previously described [70], 



with minor modifications. The leaves were ground in liquid 
nitrogen, and 30 mg of leaf tissue was homogenized with 
300 \iL of 0.1 M HCl for 5 min. The homogenate was cen- 
trifuged at 20,000 x g for 10 min, and the supernatant was 
centrifuged using a microconcentrator Ultra- 15 (Millipore, 
Billerica-USA) at 14,000 x g for 20 min. The ATP content 
was determined by using an ATP assay kit (Invitrogen) and 
a GloMax luminometer (Pr omega, Fitchburg, WI, USA). 

Measurement of mitochondrial respiration 

Measurements of mitochondrial respiration were per- 
formed using crude isolated mitochondria from six bio- 
logical replicates of WT and P07 plants. The isolation of 
crude mitochondria was performed as previously de- 
scribed [71], with minor modifications. All steps were 
performed at 4°C. Approximately 5 g of leaves were cut 
with a razor blade and ground in 20 mL of grinding buf- 
fer (0.3 M sucrose, 60 mM TES, 10 mM EDTA, 10 mM 
KH2PO4, 25 mM disodium pyrophosphate, 1 mM gly- 
cine, 1% (w/v) PVP-40, 1% (w/v) BSA, 50 mM sodium 
ascorbate, 20 mM cysteine, pH 8.0). The extract was fil- 
tered through two layers of 20 [xm nylon mesh and one 
layer of Miracloth (Millipore, Billerica, MA, USA) and 
centrifuged twice at 2500 x g for 5 min to remove the 
debris. The resulting supernatant was centrifuged at 
15,000 X g for 15 min, and the resulting pellet was gen- 
tly resuspended in 20 mL of wash buffer (0.3 M sucrose, 
10 mM TES, 2 mM EDTA, 10 mM KH2PO4, pH 7.5) 
and centrifuged at 15,000 x g for 15 min. The pellet 
containing the crude mitochondria was resuspended in 
200 [iL of wash buffer. The respiratory measurements 
were performed using an Oxygraph-2 k respirometer 
(Oroboros, Innsbruck, Austria) at 25°C in 2 mL of 
mitochondria-containing assay buffer (0.3 M sucrose, 
10 mM TES, 10 mM KCl, 2 mM MgS04, 5 mM KH2PO4, 
0.1% (w/v) BSA, 2 \iM EGTA, pH 7.5). Malate (5 mM) and 
glutamate (10 mM) were used as substrates to stimulate 
respiration. ADP (800 [iM) was added to measure ADP- 
dependent respiration. Oligomycin (5 [xM) was used as 
an FIFO ATP synthase inhibitor to measure oxygen con- 
sumption due to uncoupling activity, and the chemical 
ionophore carbonyl cyanide-4-(trifluoromethoxy)phenylhy- 
drazone (FCCP) (2 (iM) was used to determine the max- 
imum oxygen consumption. 

Transcriptome sequencing, assembly, and mapping 

Total RNA isolated from four biological repUcates of 
WT and P07 L3 leaves was used to create single-end 
RNASeq libraries using the Illumina TruSeq RNA Sam- 
ple Prep Kit (Illumina, San Diego, CA, USA) according 
to the manufacturers instructions. The libraries were 
sequenced in four lanes of an Illumina HiSeq 2000 for 
75 cycles. A total of 42.7 Gb of sequence was generated 
for the 8 libraries, with a minimum of 27.9 million 
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and a maximum of 113.2 million reads for each li- 
brary (Additional file 6: Table S5). The lUumina reads 
were filtered to remove adapters and low-quality reads 
(reads <70% bases with quality > Q20) using Adapter- 
Removal, the FASTX-Toolkit (http://hannonlabxshL 
edu/fastx_toolkit), and Perl scripts. After the filtering 
step, the reads were subjected to digital normalization 
using the diginorm software (https://github.com/ged-lab/ 
2012-paper-diginorm). The resulting 49.5 million reads 
were assembled with the Trinity [72] assembler, generating 
271,750 transcript contigs. High-quality reads (569 million 
reads) were mapped to Trinity contigs using Bowtie [73], 
and the RPKM [74] values were calculated. Because the 
N. tahacum genome is not yet completely sequenced 
and the scaffolds of N. henthamiana are estimated to 
cover only 79% of its genome, we used the complete 
genome sequence of Solanum lycopersicum [75], a close 
relative of N, tahacum, as a template. To identify protein- 
coding genes, all 271,750 contigs were used as queries in 
BLASTn and BLASTx searches against both the non- 
redundant set of 34,727 tomato-coding sequences (CDSs) 
and predicted protein sequences. The 271,750 contigs were 
compared to the Solgenomics tomato genome predicted 
protein database using BLASTx (E-value <le"^). Of the 
total contigs, 134,752 mapped to the tomato genome, 
resulting in the identification of 20,045 distinct orthologs 
for the tomato CDS models. As representative contigs, we 
chose the three best matching contigs with higher average 
RPKM values (P07 and WT) for each tomato CDS. An F- 
test was used to calculate whether the samples had equal or 
unequal variances, and the appropriate Student t-test was 
then applied to calculate the significant difference between 
the average RPKM of P07 versus WT in a contig-wise man- 
ner. The contigs identified as differentially expressed were 
selected for the next steps. The fold change value for the 
representative transcript was calculated as the ratio be- 
tween the average RPKM values in P07 and WT. The sub- 
cellular annotation of the defined contigs was determined 
by BLASTx queries against the TAIR [46] protein database 
(release 10) and TargetP [47]. Gene function data were ob- 
tained from the existing annotations of the tomato genome. 
The tomato protein dataset was annotated against UniProt 
[44] and mapped against the COG (Eukaryotic Clusters of 
Orthologous Groups) [45] database with BLASTp, and 
COG functional categories were assigned only if the two 
best hits had the same COG function. Stress-related genes 
were mapped using BLASTx against the tomato proteome 
in the Stress Responsive Transcription Factor Database 
(STIFDB) [76], and annotations were assigned to the corre- 
sponding representative tobacco contig. 

Accession numbers 

Sequence data from this article can be found in the Gen- 
Bank/EMBL database under the following accession 



numbers: A^UCPl (NM_1 15271.4), Actinl (EU938079.1), 
NADH-DeH (Y09109.1), NADH-IS2 (M77225.1), NADH- 
IS7 (L16810.1), GRPe chaperone (AF098636.1), A^MCU 
(NP_1 19075), MmMCU (XP_0065 13531), HsMCU (NP_ 
612366), A/:CCX3 (NP_566474). 

Transcriptome sequencing data are available from 
Sequence Read Archive (SRA), which is accessible through 
NCBI BioProject ID: PRJNA211804 under the experiment 
ID SUB287723. 

Additional files 



Additional file 1: Figure 51. Tobacco plants overexpressing AtUCPl 
show an increased expression of nuclear-encoded mitochondrial genes. 

Additional file 2: Table SI. Classification of mitochondrial genes 
upregulated in P07 compared to V\^. The classification was based on a 
TAIR annotation of protein subcellular prediction. 

Additional file 3: Table S2. Predicted mitochondrial genes upregulated 
with > 2-fold change in P07 compared with V\^. The prediction was 
based on TargetP. 

Additional file 4: Table S3. Classification of selected genes involved in 
different metabolic processes and significantly upregulated in P07 
compared with V\^. The classification was based on COG. 

Additional file 5: Table S4. Abiotic stress-responsive genes upregulated 
in P07 compared with WT. The classification was based on STIFDB and shows 
the most upregulated genes involved in different types of abiotic stresses. 

Additional file 6: Table S5. Summary of RNASeq data filtered and 
mapped with Bowtie. 
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